On the variational homotopy perturbation method for nonlinear oscillators 
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In this paper we discuss a recent application of a variational homotopy perturbation method to 
rather simple nonlinear oscillators . We show that the main equations are inconsistent and for that 
reason the results may be of scarce utility. 
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C$ . I. INTRODUCTION 

Q\ , There has recently been great interest in developing simple solutions to textbook models of nonlinear oscillators [l|-|4| 
i— - n , (and references therein). However, some of them are of questionable utility as shown, for example, by Rajendran et 
alQ who concluded that He's calculations of the limit cycle of the van der Pol oscillatorjl] "contain several errors 
^ pH which once rectified make the method inapplicable to it" . I have disclosed several inconsistencies in a paper by Ren 



and He Q and even proposed how to tidy up and improve their calculations Q . 

Here I discuss a recent application of a variational homotopy perturbation method to rather simple nonlinear 
oscillators In Sec. HI] I analyse their results and in Sec. Mil draw conclusions. 
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Q\ ■ II. VARIATIONAL HOMOTOPY PERTURBATION METHOD FOR NONLINEAR OSCILLATORS 

oo ; 

' Akbarzade and Langari[4| were interested in equations of the form 

o ' 

A{u)-f{r)=L{u) + N(u)-f{r) = (1) 

where L and N are the linear and nonlinear parts of the operator A and u is the solution. They proposed the 
"homotopy perturbation structure" 

H(v,p) = (1 - P )[L(v) - L(u )} + P [A(v) - f(r)] = (2) 

where p is an embedding parameter (dummy perturbation parameter in the language of the well known perturbation 
theory) and uq is the first approximation that satisfies the boundary conditions. 

They expanded the solution in p-power series v = vq + i>ip + V2P 2 + ■ . ■ and obtained the solution to Eq. (QJ as 
u = V{) + vi + v 2 + • • • provided that the series converges for p = 1 . 

In particular, the authors concentrated in nonlinear oscillators of the form 

u" + ulu + ef(u) = (3) 
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where / is a nonlinear function of u" , u' and u, and considered the "variational functional" |4J (and references therein) 

ft r 



J{u) 



X -v! 2 + -wgu 2 + eF(u) 



dt 



(4) 



where dF/du = f. Note that I have corrected a misprint in the authors' Eq. (9). Obviously, J(u) is minus the well 
known action integral Q for a particular time interval. 

In order to introduce the basic idea the authors first modified the well-known Duffing equation 



u" + u + ei? = 0, u(0) = A, u'(0) = 



(5) 



u" + lo 2 u + p [eu 3 + (1 - oj 2 )u] = 



and derived the perturbation equations of order cero 

u'q + uj 2 uq = 

and first order 

," i , ,2 



where 



satisfies the boundary conditions and 



u" + to ui + euq + (1 - uj )uq = 



Uo(t) = Acos(u)t) 



ui(0) = u[(0) = 



According the the authors "w will be identified from the variational formulation for ui, which reads" 



J(ui) = 



— -uf + -^U 2 u\ + (1 — U0 2 )u Q Ul + 6UqUi) 



2tt 

dt, T = — 

UJ 



They argued that the simplest trial function is^l\ 

u\ = B 



cos(wi) cos(5wt) 
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Surprisingly, this function satisfies one of the boundary conditions w'i(0) = but not the other one because ui(0) = 
2B/3 = leads to the unwanted trivial solution. 

From the variational conditions dJ/dB = and d J / duj = the authors obtained 



1 + -fA 2 , B = 
4 



(13) 



Although the estimated value of ui is reasonable, the result B = leads to the trivial solution u\ = that restricts 
considerably the practical utility of the approach. 

In order to improve the results the authors proposed the correction 
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_ / , . cos(3wi)\ / cos (3ut) cos(5wi)\ , 
ui = Si kos (wt) + B 3 f j (14) 

and from the variational conditions dJ/dBi = 0, dJ/dBa = 0, dJ/dui = they obtained the frequency 
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They did not show the coefficients; I obtained 



\/510237p 2 + 1416576P + 984064 - 357p - 496 (15) 



A[357p-496 (c 2 -l)] 49.4 [3p - 4 (a, 2 - l)j 
Sl " 96^ ' B3 96^ (16) 



Note that u\{t) does not satisfy one of the required boundary conditions 



- 49p - 68) 

^°) = - V 16a/ ^ 



and that the approximate solution u app (t) = Uo(t) + u\(t) exhibits the wrong amplitude u app (0) = A + ui(0) . It 
therefore seems that by means of the variational approach the authors obtained a frequency for the amplitude A and 
an approximate trajectory u app (t) with a different amplitude. 

Akbarzade and Langari[4j applied the method to other textbook nonlinear oscillators and obtained approximate 
frequencies. In all the cases they chose first-order corrections ui(t) that do not satisfy the boundary condition 
iti (0) = and obtained the trivial solution u\(t) = 0. 



III. CONCLUSIONS 



Although the combination of the homotopy perturbation method and the variational principle proposed by Ak- 
barzade and Langarij4j led to reasonably approximate frequencies, one cannot take the approach seriously because 
of its inconsistencies. First, the first-order corrections to the solutions do not satisfy one of the boundary conditions 
proposed by the authors and, second, in most of the cases the resulting corrections are trivial (that is to say, they 
vanish identically). In the case where this correction does not vanish, the approximate solution exhibits an amplitude 
that is different from the one appearing in the expression for the frequency. 
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